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Abstract 



Noise induced changes in the critical and oscillatory behavior of a Prey- 
Predator system are studied using power spectrum density and Spectral Am- 
plification Factor (SAF) analysis. In the absence of external noise, the popula- 
tion densities exhibit three kinds of asymptotic behavior, namely: Absorbing 
State, Fixed Point (FP) and an Oscillatory Regime with a well defined proper 
(natural) frequency. The addition of noise destabilizes the FP phase inducing 
a transition to a new OR. Surprisingly, it is found that when a periodic signal 
is added to the control parameter, the system responds robustly, without rele- 
vant changes in its behavior. Nevertheless, the Coherent Stochastic Resonance 
phenomenon is found only at the proper frequency. Also, a method based on 
SAF allows us to locate very accurately the transition points between the 
different regimes. 
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Oscillatory behavior is ubiquitous in all aspects of Nature. For instance, time periodic 
patterns can be found in problems involving the relation between different species (competi- 
tion, predator-prey, coexistence, etc) Oscillatory behavior can be also found in spatial 
dispersion and competition of living or chemical species and even in some physiological 
situations (cardiac and circadian rythms) 0. 

Recently, and related to the interplay between oscillations and noise, the phenomenon 
of stochastic resonance (SR) has been studied in different physical, chemical and biological 
contexts f| . In particular, SR has been found to play a relevant role in several problems in 
biology: mammalian sensory systems, increment of the tactile capacity, visual perception, 
effects of low fequency and low amplitude electromagnetic fields, etc 0. Even the problem 



of coupling among SR units has been studied [10,11 



A related, albeit slightly different phenomenon, is the so called stochastic coherent reso- 
nance (SCR) []1^1 . This phenomenon, that resembles SR, corresponds to a situation where 
the system shows noise-induced coherent oscillations [|H| without an externally applied signal 
or a discrete component in its spectrum (in this aspect it differs from SR without pariodic 
forcing as discussed in [p|]). Here we analyze the effect of fluctuations on a lattice gas 
model for a prey-predator system with smart pursuit and escape. Our aim is to analyze the 
possibility that a SCR-like phenomenon can occur in such a system. 

Our lattice gas model, which is a variant of the cellular automata proposed by Boccara 



et. al ||15|| , is defined as follows: a lattice site can be either empty or occupied either by 
a prey or a predator. Double occupancy of lattice sites is forbidden. The system evolves 
according to consecutive cycles: i) competition of species and ii) escape-pursuit dynamics. 

The competition rules are as follows: (ajPreys have an offspring occupying an empty 
next neighbor site with probability Bpu (Birth Probability of Preys) in case of absence of 
predators within their Vrh (visual range of preys), (b) Predators can eat a prey that is 
in their Mrp (Movement Range of Predators) with probability Dph (Death Probability of 
Preys), (c) Predators who have already eaten a prey can produce an offspring in the site 
occupied previously by the eaten prey, with probability Bpp (Birth Probability of Predators). 



(d) Predators can suddenly die with probability Dpp (Death Probability of Predators). 

The rules for an escape-pursuit process are as follows: (a) Preys calculate the gradient of 
the density of predators in their Vrh and move into an empty site in the opposite direction. 
(b) Predators calculate the gradient of the density of preys in their Vrp (Visual Range of 
Predators) and move into an empty site in that direction. 

We have restricted ourselves to investigating the dependence of the system on the preda- 
tor birth probability (Bpp), while the remaining parameters are kept constant, namely: 
Mpp = Vrp = Vrh = 1, -Dpp = 0.25, Dpp = 1. and Bpp = 0.5. The model is studied 
within a mean field approach. Further details on the model and results of Monte Carlo 
simulations have already been published . 



With the above considerations, the mean field equations for the system can be derived 
evaluating the rates of all processes which may change the species densities. Then, it is 
obtained 

d tPP = p P -(B-D) (1) 
d tPH = PH ■ {A - (B + C) ■ [l - (1 - PP f M ^ 2 - 1 ] } (2) 

where pp and pp are the global predator and prey densities, respectively. Also the rates are 
given by: (i) The rate of prey reproduction A = B PH ■ (1 — Z)p P )( 2V «- ff+1 ) 2 ~ :L [l — (p P + Ph) 8 ] 
which comprises the probability of a prey to have an offspring (Bpp:) in a neighboring empty 
state (first bracket) provided the absence of predators in the neighborhood of the progenitor 
(second bracket), (ii) The rate of predator reproduction B = (1 — D PP ) ■ B PP • D PH ■ [1 — 
(1 — PhY 2Mrp+1 ^~ 1 ], which comprises the probability of the predator to survive (1 — Dpp), 
the birth probability of predator (B PP ), the death probability of a prey (D PH ) and the 
probability of existence of a prey in the visual range of a predator (last term), (iii) The 
probability of a predator to catch a prey C — (1 — Dpp) ■ (1 — Bpp) ■ Dpp • [1 — (1 — 
equal to the previous term B, except for the fact that the predator does 
not have an offspring and consequently Bpp has been replaced by (1 — Bpp). (iv) Finally, 
the dying probability of predators is D = Dpp. 



In order to study the influence of external perturbations on the system, it is assumed 
that the control parameter, Bpp, is time dependent and has the form 

Bpp(t) = Bp Po + QZ(t)+ecoa(n e t), (3) 

where is a normalized Gaussian white noise of intensity Q and e is the amplitude of a 
periodic external signal with frequency Q £ . In any biological system a parameter such as the 
Birth Probability of Predators (or any other) will not be constant in time. Effectively, the 
environment competition, climate variability, etc, will alter the characteristics of the birth 
rate of all species. Although it is true that all the parameters of this system will change in 
time, as a first step in our study, we simplify the analysis by only considering modifications 
on the Bpp parameter. 

Clearly, the choice of a Gaussian noise implies the possibility of having negative values of 
Bpp(t). Such a case has no physical meaning but, as the width of the Gaussian distribution 
is too small, the small negative tail does not affect our results. 

Initially, the behavior of the system has been studied in the absence of any external 
perturbation, i.e. fixed B PP . It was found that the densities of species, pp and pu show the 
following asymptotic behavior depending on the value of Bpp: (i) If Bpp < 1/3 the system 
evolves toward the extinction of predators, namely pp = 0, and p# = 1. Such a state is an 
"Absorbing Regime" (AR), because the system can not escape from it. (ii) If Bpp > 1/3 
the final state of the system is a steady regime with coexistence of prey and predators. 
However within this regime, dynamic behavior changes according the value of Bpp. Thus, 
if Bpp is near the AR phase, pp and pu reach constant values, say a "Fixed Point" (FP) 
regime (see Fig. |]). On the other hand, for high enough values of Bpp, the system enters an 
"Oscillatory Regime" (OR), since both populations exhibit self-sustained oscillations (Fig. 
[1]). It should be noted that the proper frequency of the system is independent of B PP . 

Hence, Bpp is the control parameter governing the system dynamics, while pp and pn 
are order parameters. Then, it is possible to infer that there are critical values Bp c p and 
Bp c P separating both the AR and the FP, and FP and OR, respectively. 
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The discrete versions of Eqs. ([!]) and @ have been solved in order to: (i) determine if 
the system shows the SCR phenomenon; (ii) provide methods for the accurate localization 
of the critical values. 

Time series of population densities, puif) and p#(i), were computed recursively from 
the discrete mean field equations. Then, the respective temporal self-correlation functions 
Kh{t) and Kp{r), given by 



were calculated. Here ()t represents the time average. It is well known that the power 
spectrum density (psd) is just the Fourier Transform of such functions Sh(u) = ^[Kuij)] 
and Sp(uS) = jF[i^ P (r)], respectively. Then, the SNR is readily obtained as 



where Sfj (Sp ) is the output power at the frequency of the signal considered, and (Sp) 
is the output power of the noisy background at the same frequency, both of them obtained 
from the prey (predator) time-series, respectively. From now on, we will present the results 
for the predator time-series only, as the prey ones shows qualitatively the same behavior. 
Also, the subscripts will be omitted to simplify the notation. 

It is worthwhile to remark here, that the analysis was made considering the proper 
frequency of the system, Q* = 0.03. At this frequency R reaches a significant value, in 
contrast to its value at the forcing frequency (Q £ ) where R is negligible. We also verified 
that the periodic term in Eq. ([3]), does not introduce changes in the dynamics of the 
system. Thus, we will not consider in this work any oscillating perturbations on the control 
parameter. 

Figures |2| and |3| show the psd when the system is in the FP and OR regimes, respectively. 
Qualitatively, the results can be understood as follows (see Fig. ^]): if the system is in 
FP, the psd at the proper frequency is zero when Q = 0. For very low noise strengths 
(Q < 10~ 4 ), S(Q*) is an increasing function of Q. However, the situation changes when the 



K H (r) = (p H (t)p H (t + r)) 



K P {r) = (p P (t)p P (t + T)) t . 



(4) 
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noise is strong enough (Q ~ 10~ 3 ). Here, the peak at the proper frequency is hidden by the 
noisy background. Then, it is clear that the function R(Q) would have a maximum in the 
considered range of Q. In contrast, within the OR (see Fig. ^) S(Q*) has a finite (non-zero) 
value in the limit Q — > 0. As the noise strength grows, such a peak decreases and broadens. 
Then, for this case the SNR of the system as a function of noise strength is expected to be 
a monotonic decreasing function of such a parameter. 

In the following, a method to locate the phase transition between FP and OR regimes, 
by means of the calculation of the SNR, is developed. Figure f| shows plots of the SNR as 
a function of Q for different values of Bpp Q . The SCR phenomenon is only apparent for 
values of B PPo < 0.47. Such a dependence on Q is the main characteristic of the system 
for the FP regime. When B PPo > 0.47 the SNR is a decreasing function of Q, which is the 
expected behavior when the system is in OR. The value of the parameter Bp c P where the 
SCR phenomenon disappears can readily be identified with the critical point. Hence, our 
best estimate for such a transition point is given by B 2 P C P = 0.47 ± 0.01. 

Figure ^| shows the dependence of SNR as a function of B PPo , for some fixed values 
of Q. In the absence of noise the critical point can easily be obtained from the plot (Fig. 
|S](a)), yielding Bp P = 0.47 ±0.01, in agreement with the value obtained through the method 
discussed above. 

A surprising consequence of the addition of noise is that, as can be seen in Figs. H(b)-(e), 
even for control parameter values B PPo < 0.47 the system oscillates at the proper frequency. 
This indicates that the FP phase is unstable, and even a small perturbation drives it into 
an oscillatory behavior. Also, the plots of R vs B PPo exhibit linear behavior for lower values 
of the parameter, e.g B PPo < 0.46. Extrapolations by linear regression to the limit R — > 
give Bp c p a = 0.330 ± 0.005 which is the critical point for the irreversible transition from 
the coexistence regime to the AS. This figure is in excellent agreement with our previous 
estimate B P C P = 1/3 obtained simply solving the mean field equations in the absence of 
noise, as discussed earlier. 

Subsequently, for small noise values Q < 10~ 4 (Figs. ||(b)-(c)), the onset of a shoulder 
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is clearly observed, say for Bpp a > 0.47. As will be discussed later, this can be associated 
to another phase transition between intrinsically different oscillatory regimes. 

In the search of a method to characterize these two different oscillating phases we make 
use of the fact that, although this is an extremely non-linear system, the time series pn and 
p P are composed only by a single frequency (the proper one) and its harmonics. Supported 
by this numerical evidence, it is assumed that for long times, the solution of Eqs. ([!]) and 
(@) approaches a periodic function, which will be of Floquet-type. Then, the corresponding 
stationary self-correlation function (K as (r)), can be decomposed in a Fourier series, i.e. 

oo oo 

K(t) = K sin(nfir) + £ B n cos(nfir). (6) 

n=l n=0 

It is expected that in the FP regime, all the coefficients (with the exception of -Bo) will vanish. 
Also, in OR, at least the terms A\ and B\ (related with the amplitude of oscillations) must 
be non zero. Hence, measuring the value of those coefficients as a function of Q and Bpp it 
may be possible to find the phase transition. 

Once again, this is a SCR phenomenon, as analyzed within the framework used in Ref. 
[T7|] , where the characterization of SR is done by means of the Spectral Amplification Factor 
(SAF). Following this course, we define W, the degree of oscillation (this will indicate how 
much inside the OR is the system), as the inner product between K as (r) and the first 
harmonic of Eq. (§) , which gives 

U = A\ + B\. (7) 

Figure |6] shows plots of W as a function of the control parameter Bpp a , in the absence 
(filled circles) and presence (empty circles) of noise. It is easy to distinguish the following: 
(i) In the absence of noise the FP regime is clearly observed in the Figure, and a value Bp c p = 
0.47 is found, in excellent agreement with the results previously discussed in connection with 
the behavior of SNR vs. Bpp Q (Fig. |5|). (ii) In the presence of noise it can be seen that the 
value of W increases for Bpp a < 0.45. Then, a new oscillatory regime is apparent where W 
is independent of Bpp a , and such a behavior arises instead of the FP phase, because of the 
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applied low-intensity noise. This phenomenon can be understood in terms of a noise induced 
phase transition. On the other hand, the usual OR where W monotonically increases with 
Bpp , is obtained for Bpp Q > 0.45. In fact, the value of Bpp Q for which the system enters 



into the usual OR, suffers a shift due to the presence of noise fLo] . Effectively, it changes 
from B PPo (Q = 0) = 0.47 to B PPo (Q = lO" 5 ) = 0.45. 

The former study of W behavior allows us to give further support to our discussion of the 
results shown in Fig. ||. Substracting the linear behavior, the following estimates for the 
phase transition are obtained: B PPo (Q = 10~ 5 ) = 0.461 ± 0.005, and B PPo (Q = 10~ 4 ) = 
0.453 ± 0.005. Finally, when Q increases (Figs. ||(d) and|^(e)), the linear dependence on 
Bpp holds for a wider range in such a variable and consequently, the phase transition is 
more difficult to detect. This is because the system is strongly oscillatory even if (in absence 
of noise) Bpp a is within the FP regime. In fact, the linear tendency in the function R(Bpp ) 
observed for small values of Q is broken when the parameter is set within the OR. Also, 
when the OR is masked by the response to a larger noise, the linear dependence is not broken 
but the critical point can not be located anymore. Furthermore, from Figs. |5|(a) and ^| one 
observes that in the absence of noise the transition FP <-> OR is abrupt resembling first 
order behavior. However, when the noise is added the transition between the two oscillatory 
regimes not only becomes shifted, but also rounded (Figs. |(b) and ||). 

Summarizing, the effect caused by external noise on a Prey-Predator system, with smart 
pursuit evasion is studied by means of SNR and Degree of Oscillation techniques. In order to 
perform the study, the behavior of the system is analyzed when a small amount of white noise 
is applied over a control parameter (Bpp ). It is found that if (in the absence of noise) such 
a parameter is in FP, a Coherent Stochastic Resonance phenomenon arises. In fact, only the 
proper frequency becomes strongly enhanced, while any externally applied oscillatory signal 
is not amplified. Nevertheless, such an enhancement arises only for a (intermediate) range 
of noise strength (Q ~ 10 -6 ) and it vanishes in two distinct limiting cases: (i) when noise 
is not applied and (ii) when the noise is too large (Q > 10~ 3 ). It should be remarked that 
this resonant phenomenon disappears if the system is driven to the OR, tuning the control 
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parameter. Then, it is possible to locate, looking at the dependence of SNR as a function of 
noise, the critical point at which the phase transition between the FP and the OR occurs. 

In the presence of noise, the FP phase is replaced by a new oscillatory regime where the 
amplitude of oscillations remains constant, allowing us to assert that the FP phase is less 
stable than the OR one. The dependence of SNR as a function of B PPo shows linear behavior 
and the extrapolation of such a line intersects the SNR axis at Bpp Q = 0.330 ± 0.005. This 
value is the critical point where the irreversible phase transition to the AS occurs. Thus, 
this method allows us to find both phase transitions: (i) the irreversible one between th 
Coexistence Regime and the AS, and in the absence of noise (ii) the reversible one between 
the FP and the OR regimes, respectively. 

We have also defined the Degree of Oscillation as the corresponding terms of the first 
harmonic in the Floquet-type asymptotic expansion of the self-correlation function of the 
population time series. This characterization, provides us with an alternative method to 
locate the phase transition between FP and OR. 

Furthermore, the abupt (first-order like) transition between the FP and the OR in ab- 
sence of noise becomes clearly rounded (second-order like) when noise is applied. These 
results point out that the effect of noise is not merely restricted to the shift of the critical 
point, but instead that the whole nature of the transition changes. 

We expect that the findings reported in this work will contribute to the understanding 
of resonant effects and critical behavior in actual competitive population systems, since in 
nature they are exposed to different sources of external noise. 
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FIGURES 




FIG. 1. Phase diagram of the prey-predator system as obtained solving the mean- field equations 
in absence of noise, (a) Plot of Dpp versus Bpp showing the critical curve (solid line) for the 
irreversible phase transition between the prey-predator coexistence regime and the absorbing state 
where predators become extinct. The dotted line shows the critical curve for the transition between 
the fixed point (FP) regime and the oscillatory behavior (OB). Figures (b) and (c) show the 
dependence of the population densities for the different regimes of the system, namely the FP and 
the OB, respectively. 
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FIG. 2. Plots showing the dependence of the psd when the noise is modified. In absence of 
noise the system is in the Fixed Point regime since Bpp Q = 0.42 is kept constant. In each graph, 
we varied Q, the values being: (a) Q = l(T 8 ,(b) Q = 1(T 2 , (c) Q = 4.45 • 10" 2 , (d) Q = 10 -1 . 
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FIG. 3. Plots of 5 versus uj corresponding to the time series of the predator density. In this 
case, when Q = the system is in OR. The results are for: (a) Q = 0, (b) Q = 1(T 4 , (c) Q = 1(T 2 
and (d) Q = 1CT 1 , keeping Bpp D = 0.70 constant. 
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FIG. 4. Plots of the Output SNR as a function of noise strength, Q, obtained for different values 



log, "1(0 



of Bpp Q . A SCR phenomenon is apparent when Bpp Q is within FP regime, but it disappears if the 
system is within OR. (a) B PPo = 0.42, (b) Q = B PPo = 0.45, (f) B PPo = 0.47, (d) B PPo = 0.50. 
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FIG. 5. Plots of R as a function of the control parameter Bpp Q . In each plot, the noise strength 
is kept constant, as follows: (a) Q = 0, (b) Q = 1(T 5 , (c) Q = 1(T 4 , (d) Q = 1(T 3 and (e) 
Q = 10- 2 . 
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FIG. 6. Plots of the Degree of Oscillation as a function of Bpp Q , obtained for Q = (a), and 
Q = 10~ 5 (b). In both cases, the phase transition between FP and OR can be observed. Notice 
that such a transition is shifted in the presence of a small noise strength. 
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